arXiv:l508.03114v2 [cond-mat.mes-hall] 29Jul2016 


Chiral anomaly in Dirac semimetals due to dislocations 
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The dislocation in Dirac semimetal carries an emergent magnetic flux parallel to the dislocation 
axis. We show that due to the emergent magnetic held the dislocation accommodates a single fermion 
massless mode of the corresponding low-energy one-particle Hamiltonian. The mode is propagating 
along the dislocation with its spin directed parallel to the dislocation axis. In agreement with the 
chiral anomaly observed in Dirac semimetals, an external electric held results in the spectral how 
of the one-particle Hamiltonian, in pumping of the fermionic quasiparticles out from vacuum, and 
in creating a nonzero axial (chiral) charge in the vicinity of the dislocation. 

PACS numbers: 75.47.-m,03.65.Vf,73.43.-f 


I. INTRODUCTION AND MOTIVATION 

The Dirac semimetals are novel materials that have been discovered recently (NaaBi and CdaAs 2 [IH2I). A possible 
appearance of Dirac semimetals in the other systems (for example, ZrTe 5 01 [5], and BDSea 0 ) was also discussed. 
In Dirac semimetals the fermionic quasiparticles propagate according to the low energy action that has an emergent 
relativistic symmetry. Both in NasBi and CdsAs 2 there exist two Fermi points ±K^°\ At each Fermi point the pair 
of left-handed and right-handed fermions appears. The Dirac semimetals represent an arena for the observation of 
various effects specific for the high energy physics. In particular, the effects of chiral anomaly play an important role 
in physics of these materials 0EHHI- 

In the Weyl semimetals, which were also discovered recently (in particular, TaAs [12]) one of the two Fermi points 
hosts a right-handed Weyl fermion while another Fermi point hosts a left-handed Weyl fermion. Various relativistic 
effects were discussed in Weyl and Dirac semimetals already before their experimental discoverv |13U25j . 

In 0] the experimental observation of the chiral anomaly in ZrTes was reported as measured through their contri¬ 
butions to the conductance of the sample. It has been shown, that in the presence of parallel external magnetic field 
and external electric field the chiral anomaly leads to the appearance of nonzero chiral density and, correspondingly, 
a nonzero chiral chemical potential. This work was followed be a number of papers, where the experimental detection 
of chiral anomaly was reported in different Dirac and Weyl semimetals (see [55] and references therein). 

Similarly to graphene [27U34] the fermionic quasiparticles in Dirac and Weyl semimetals experience emergent gauge 
field and emergent gravity in the presence of elastic deformations of the atomic lattices (see, for example, [55M58] and 
references therein). In this paper we will concentrate on dislocations in the crystalline order of the atomic lattice, 
which are particularly interesting cases of the elastic deformations of the ion crystal lattice [39l 140] . The dislocation 
is a line-like defect characterized by the Burgers vector b which determines the physical displacement of the atomic 
lattices along the dislocation. The vector b is a global characteristic of the dislocation because it is a constant quantity 
over the entire length of the dislocation. In rough terms, one may imagine the dislocation as a vortex which possesses 
a fixed “vorticity” given by the Burgers vector b. The extreme examples of the dislocations are the screw dislocation 
(shown in Fig. |T|) and the edge dislocation (illustrated in Fig. [2]) for which the corresponding Burgers vectors are 
parallel and, respectively, perpendicular to dislocations’ axes n. There are other types of the dislocations lying in 
between these two extreme cases. 

In [38j the effect of the dislocation on the geometry experienced by fermionic quasiparticles in Dirac semimetals 
was considered for the first time. Aharonov-Bohm effect and Stodolsky effect (the latter effect describes a correction 
to the Aharonov-Bohm effect due to torsion) were investigated for the scattering of the quasiparticles on dislocations. 
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FIG. 1. Illustration of the screw dislocation of the atomic lattice with the Burgers vector b parallel to the axis n of the 
dislocation (the green line). The semitransparent plane points out to the region where the atomic planes experience a shift. 



FIG. 2. Illustration of the edge dislocation with the Burgers vector b perpendicular to the dislocation axis n (the blue line). 
The semitransparent plane shows the extra half-plane of ions introduced in the crystal. 


Besides, basing on an obvious analogy with the results of [H] it was proposed, that the dislocation (that carries an 
emergent magnetic flux) becomes the source of chiral anomaly and chiral magnetic effect. This occurs because the 
dislocation carries emergent magnetic field. Therefore, it was argued, that the chiral anomaly and chiral magnetic 
effect occur without any external magnetic field. According to [38] the contribution of topology to magnetic flux <I> is 
equal to the scalar product K^b, where b is the Burgers vector. There may also appear the contribution to the flux 
$ proportional to the tensor of elastic deformations caused by the dislocation with the coefficients of proportionality 
that are analogous to the Gruneisen parameter of graphene. The emergent magnetic flux is associated with emergent 
magnetic held. 

H\x) = <S> J dy i S ( - 3 \x — y), (1) 

where the integral is taken along the dislocation. The appearance of the delta-function in Eq. ([!]) in the low- 
temperature theory corresponds to the fact that the emergent magnetic flux is localized within the dislocation core 
of radius £ ~ |b|, where b is the Burgers vector of the dislocation. In [38] the simple model of the dislocation was 
used, in which it is represented as a tube of size £ with the emergent magnetic held inside it. 

The further examination of the mentioned above problem has led us to the conclusion, that the naive application 
of the pattern of chiral anomaly discussed in [4] to the case, when the magnetic held is emergent and is caused 
by dislocations, has certain restrictions. Strictly speaking, the mentioned above model of the fermionic excitations 
and chiral anomaly within the dislocation may be applied to the investigation of real materials only if the emergent 
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magnetic flux of the low energy field theory is distributed within the area of size £ essentially larger, than the 
interatomic distance a while the emergent magnetic flux of the dislocation is essentially larger, than 2i r. In this 
situation we are formally able to use the low energy field theory for the description of fermionic excitations inside the 
dislocation core. This occurs for the strong dislocations with sufficiently large values of the Burgers vector b, when 
the crystal lattice is distorted considerably (or in the case, when many parallel dislocations with small values of the 
Burgers vectors are located close to each other). In this case the dislocation core size £ ~ |b| a is much larger than 
the interatomic distance a. 

In the present paper we consider the opposite situation, when in Dirac semimetals the values of Burgers vector 
are relatively small, so that the magnetic flux at the dislocation is smaller than 2n or around 2n. In this situation 
the crystal structure is not violated strongly, so that the dislocation core size is, presumably, of the order of the 
interatomic distance £ ~ a. The low energy theory is developed for the states with the typical values of momenta 
much smaller, that 1/a. Therefore, in this case the states localized within the dislocation core cannot be described 
by the field theory. In order to describe such states the microscopic theory is to be applied. 

It appears, that the anomaly in the right- and left-handed quasiparticle currents is given by 

(d,f R>L )=±^EB or ( d ^) = ^EB, (2) 

where the upper and lower signs in the first equation correspond to the right-handed and the left-handed quasiparticles, 
respectively. The important feature of Eq. (|2j) is that the effective magnetic field B - contributing to the anomaly ([2]) 
differs from the emergent magnetic field at the dislocation H as given in Eq. (jTJ) . 

The basic reason for the difference between the emergent magnetic field H and the effective magnetic field B is 
that the emergent magnetic field H of the dislocation has a small (of the order of unity or even smaller) magnetic 
flux <E>. In this case the contribution to both mentioned effects is given by the single fermionic mode (related to the 
zero mode of the one-particle Hamiltonian) propagating along the dislocation rather than by a large ensemble of the 
lowest Landau modes with a huge degeneracy factor. To be more precise, the effective magnetic field B is expressed 
through the probability density corresponding to the zero mode of the one-particle Hamiltonian in the background 
of the emergent magnetic field H due to the dislocation. The appearance of the propagating (zero) mode at the 
dislocation is a natural effect, which is known to exist in topological insulators with lattice dislocations [45] . 

The contribution of the individual zero mode to the chiral anomaly © can be described with the help of the effective 
magnetic field B, which carries a unit elementary magnetic flux contrary to the original emergent magnetic field if, 
which may have an arbitrary (but still small) value of the total flux <I>. The effective field B is localized in the wide 
area of linear size where l/£ 0 is the infrared cutoff of the considered field theoretical low energy approximation 
(below we argue that £o may be identified with the mean free path of the quasiparticles. For example, in CdaAs 2 the 
mean free path is £o ~ 200 /jm 0). 

In this paper we demonstrate that for the straight screw dislocation directed, for example, along the symmetry axis 
of the crystal the emergent magnetic flux associated with the emergent field H is given by 

$=(ifW.6) + ^(n-6), (3) 

where the first term is of the topological origin ®Sj while the second term is not topological (here (3 is an analogue 
of the Gruneisen parameter of graphene [35]). The magnetic field associated with the flux ([3| is localized within the 
dislocation core of a typical size £ ~ a, where a ~ Iran is a typical interatomic distance. In Eq. ^ the vector K ^ 
encodes positions of the Fermi points k = in momentum space, and n is the direction of the dislocation axis. 

For the straight screw dislocation vectors , b and n in Eq. (|3| are parallel to each other. 

We will discuss effects, which appear due to the interplay between quantum anomaly and dislocations in the crystal 
structure of Dirac semimetals. A fermion excitation is affected by the dislocation, in particular, via the mentioned 
above intrinsic magnetic field, which is localized in a spacial vicinity of the dislocation and is directed along the axis 
of the dislocation. In principle, the emergent magnetic fields corresponding to different Weyl fermions (that belong to 
different Fermi points and/or have different chiralities) differ from each other. However, there exists an approximation, 
in which those emergent fields H have the same absolute values, but opposite directions for the two Fermi points 
If this approximation is not violated strongly (which is the general case) the signs of the emergent magnetic 
fluxes experienced by the quasiparticles living near to the Fermi points ±K i '°' 1 are opposite. In the Dirac semimetal 
both right- and left-handed fermion excitations are present in each (of the two) Dirac cone, therefore in this case we 
have a standard effective magnetic field B(x) acting on the right- and left-handed fermions at one Dirac cone and the 
magnetic field —B(x) acting on the right- and left-handed fermions at another cone. These fields enter expression for 
the anomaly Eq. ©• 

If now one applies an external static electric field E along the axis of the dislocation, then the quantum anomaly 
will generate the chiral charge at a rate proportional to the scalar product EB. The generated chiral charge will 
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dissipate, both due to the chiral-changing processes inside the region of size £o around the dislocation and due to 
the spatial diffusion of the chiral charge. Next we notice that the equilibrium distribution of the chiral charge - 
which can effectively be described by a spatially nonconstant but otherwise static chiral chemical potential - is 
subjected to the intrinsic magnetic field of the dislocation itself. The chemical potential (j,§ is distributed around 
the dislocation with the characteristic length Ly (for example, in CdaAs 2 this length is of the order of Ly ~ 2 /im.) 
The chirally imbalanced matter in the presence of magnetic field generates dissipationless electric current directed 
along the dislocation and concentrated in the spatial vicinity around it. Therefore, the intrinsic magnetic field of 
the dislocation would lead to a spatially-dependent (negative) magnetoresistance around the dislocation. Similar 
arguments were used in Ref. [5] to experimentally investigate effects of the chiral anomaly in ZrTes in the presence of 
external magnetic field. 

The paper is organized as follows. In Sect. [Tl]we recall briefly general theory of quasi-relativistic fermions in Dirac 
semimetals in the presence of elastic deformations which leads both to the emergent gauge field and to the emergent 
gravity (the latter is described by an emergent vielbein ; 3Tj). In Section III we discuss these effects focusing on 
dislocations, partially following Ref. [3S]. In Sect. IV we consider the zero modes of the one-particle Hamiltonian and 
demonstrate, that there always exists a single mode with the definite spin directed along the emergent magnetic flux, 
which is localized in a wide area around the dislocation. In Sect. [V]we show that the spectral flow along the branch 
of spectrum (that crosses zero at the mentioned zero mode) gives rise to the anomalies in quasiparticle currents: in a 
Dirac semimetal the chiral anomaly appears. For the sake of simplicity, these results are discussed first for the strait 
dislocation, is directed along the symmetry axis z of the crystal that coincides with the direction of the Fermi p oint 
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in momentum space. We extend our results to the case of strait dislocations with arbitrary direction in Sect. 
Then in Sect. VII we discuss the generation of the chiral charge via the chiral anomaly ([2]) due to the interplay between 
an external electric field and the internal magnetic field of the dislocation. The last section is devoted to discussions 
and to our conclusions. 


II. RELATIVISTIC FERMIONS IN DIRAC SEMIMETAL 


The Dirac semimetal possesses two cones, each of which hosts one right-handed and one left-handed Weyl fermion. 
In the presence of elastic deformations caused by the dislocation the action for a right-handed and left-handed Weyl 
fermions near a given Fermi point are, respectively, as follows [35]: 


J *x\e\[*u%( X )o>W - [D^]ie ^{x)a b ^] : 

(4) 

Sl=^ J d 4 x\e\[&ieg(x)a b V M & - [D^]ie^(x)a b ^], 

(5) 

where 


iVf, = + A m (: r) 

(6) 

is the covariant derivative corresponding to the emergent t/(l) gauge held A^, a 0 = a° = 
a = 1,2,3 are the Pauli matrices. The currents of the right- and left-handed quasiparticles are, 

1, and (7° = — a a with 
respectively, as follows: 

j r = ^ie%(x)a b ' F, 

= ^ie^{x)a b \& . 

(7) 

(8) 

Throughout this paper the internal ^0(3, 1) indices are denoted by Latin letters o,6, c,... from the beginning of the 
alphabet while the space-time indices are denoted by Greek letters or by Latin letters i,j, k, ... from the middle of the 
alphabet. 

The vierbein held = e%(x) is a 4 x 4 matrix which carries all essential information about anisotropy and the 


elastic deformations (caused, for example, by a dislocation) of the ion lattice of the Dirac crystal. It is convenient to 
introduce the inverse of the inverse vierbein held, e“ = e“(a;), defined, naturally, as follows: 


e a( x ) e v( x ) = ■ (9) 

In our paper we always assume that the deformations are small so that the determinant of the vierbein held 

|e| = det(e“), 


never vanishes. 


( 10 ) 
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In the absence of elastic deformations the fields entering the actions Q and ([5]) are simplified. In this case the 
emergent gauge field A M vanishes. 

In the absence of elastic deformations the vierbein can be chosen in a diagonal form, 


b (o).m 


( v 

0 
0 

V o 


-1 

F 


0 0 0 \ 

V - 1 / 3 0 0 

0 v~ x ! 3 0 

0 0 V 2 / 3 ) 


( 11 ) 


where the parameter v ^ 1 reflects the fact that the experimentally studied Dirac semimetals are anisotropic materi¬ 
als DM33- It is also convenient to introduce the spatial component of the undeformed vierbein (|11[): 


=(o ),» = fi — 

-'a — J a 


V - 1 / 3 0 0 

0 v~ x l 3 0 

0 0 v 2 ! 3 


( 12 ) 


with i,a = 1,2,3. The quantity = Vpf-' with fixed i = 1,2,3 has a meaning of the anisotropic Fermi velocity in 
i-th direction. The determinant (10) in the unxleformed case is |e(°)| = vf- 

\K (0) \ 


The low-energy effective field theory (|4j), (|5| has the natural ultraviolet cutoff A uv ~ \K { j | associated with the 
positions of the Dirac cones in the momentum space. In order to determine a natural infrared cutoff we notice that 
in our field-theoretical approximation the massless quasiparticles do not interact with each other since the effective 
actions Q and (|5j) contain only bilinear fermionic terms while the gauge field A^ is a classical non-propagating field. 
Therefore, the natural infrared cutoff for our approach is A/^ = l/£ 0 , where the length £ 0 may be identified with 
the mean free path of the massless quasiparticles. Indeed at the distances of the order of the mean free path £o we 
cannot neglect interactions between the quasiparticles and their scattering on the defects of the atomic lattice which, 
in general, cannot be captured by Eqs. Q and ([5]). 

As an example, we mention that for the Dirac material Cd 3 As 2 the mean free path £ 0 was estimated in Ref. [ 8 j to be 
of the order of 200 /im. In the above formulation of the low-energy theory, the Dirac point corresponds to zero energy. 
In real situation the crystals of semimetal may have nonzero Fermi energy at the level crossing points. In particular, 
in [42) the values of Fermi energy of the order of 10 meV were reported for NasBi. In the following applications we 
assume, that in the real systems the value of Fermi energy may be neglected, or that the sample is doped in such a 
way, that the doping-induced chemical potential shifts the level crossing to the vanishing energy. 

In the upcoming sections for simplicity we restrict ourselves to the case, when the dislocation is an infinite straight 
line directed along the symmetry axis z of the crystal, which coincides with the direction of the Dirac point in 


momentum space. We will return to the more general case of an arbitrarily aligned straight dislocation in Sect. VI 


Now let us consider the case when the atomic lattice of a Dirac semimetal is elastically deformed. The deformation 
is described by the displacement vector u l which gives the displacements of the ions with respect to their positions 
with respect to the unperturbed semimetal. In the approximation of isotropic elasticity for a straight dislocation 
directed around the z = x 3 axis the displacement vector u a is given by: 


b a 

u ° = 


(13) 


where 6 is the polar angle in the plane orthogonal to the dislocation and b a is the Burgers vector. The first term in the 
right hand side of Eq. (13) is discontinuous vector function as it has a jump by b a at 8 = 0. The second, continuous 


part of displacement is given by [44] 


t (*-!-) = - 


b l 1 — 2 cr 
4-7T 1 — a 


_3 kl 


log 


\x ± \e^ e 3il x\x k 




2 R 


1-2(7 


(14) 


where x±_ = ( 24 , 2 : 2 ) are the transverse coordinates in the laboratory reference frame, x j_ are respective unit angles 
in the transverse plane and a is the Poisson ratio which is defined as the negative ratio of transverse to axial strain 
of the atomic crystal. Throughout this paper we shall work in the laboratory reference frame in which the positions 
of ions are their real 3d coordinates. 

Notice, that for the screw dislocation when the Burgers vector directed along the dislocation axis, b = (0, 0, b z ), 
the continuous part of the displacement vector vanishes, z((? ont = 0. It is worth mentioning, that while the values of 
u k ont may be large, its derivatives are small for sufficiently small b because after the differentiation the expression in 


Eq. (14) tends to zero at |ajj_| —>• 00 . 

In the presence of elastic deformations, in principle, the emergent vielbeins (as well as the emergent gauge fields) 
may differ for the left-handed and the right-handed fermions incident at the given Dirac point. 
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Let us introduce tensor of elastic deformations [32] 


u ij = <)' vP + HV, 


(15) 


where we have neglected the part quadratic in u l by assuming that the deformations are small. In general, the 
emergent vielbein around the dislocation may be expressed, up to the terms linear in displacement vector, as follows 
(see Ref. [158: for the details of the derivation): 

e a = /i(l + \lL 3 u n >) + tidku 1 t'& jh v? k 

4 = -—e°=0 

Vp J 

eo = — (1 + 

Vp 6 J 

|e| = v F (l - dm 1 - ^7 kijU lJ ) 

a,i,j,k,n= 1 , 2,3 ( 16 ) 

The emergent gauge field is given by 

A i n-V i {u-KW)) + -P ijk u,i\ ( 17 ) 

a J 

A o = ^Pojk u3k ’ = 1 , 2,3 


The tensors /3 and 7 , which are the analogues to the Gruneisen parameters in graphene, may, in principle be different 
for the right-handed and the left-handed fermions. The analogy to graphene prompts that their values could be of the 
order of unity. Notice, that in graphene the emergent electric potential Ao does not arise outside of the dislocation 
core [5? . In the same way we assume, that in the semimetal the parameters /3 0 jfc may be neglected. The reason for 
this is that the combination K + A appears as the value of momentum P , at which the one-particle Hamiltonian 
H(x, P) vanishes (one substitutes K (0 ^ + A instead of the momentum operator P): 

7£(ai, JC (0) + A(x)) =0 (18) 


As a result we expand the Hamiltonian near the floating Fermi point K ^ + A(x): 


H(x, P) = |e(*)| e k a (x)a a o P k - (k[ 0) + A k (x)) 
+Ao(a:), 


(19) 


where by the symbol o we denote the symmetric product 


AoB= -(AB + BA). 


( 20 ) 


The only possible source of Aq(x) is the noncommutativity of momentum P and coordinates. This means, that 
unlike A k with k = 1,2,3 the emergent electric potential Ao is proportional to the derivatives of the parameters 
entering H(x,P). The field A k with k = 1,2,3 is proportional to 1/a times the combination of the dimensionless 
parameters while Aq is proportional to their derivatives but it does not contain the factor 1/a. For slow varying elastic 
deformations this means that A 0 may be neglected. This consideration does not work, however, within the dislocation 
core, where physics is much more complicated. The influence of this unknown physics on the quasiparticles with small 
values of momenta (described by the action of the form of Eqs. & ©) may be taken into account through the same 
emergent fields = 0,1,2,3 and e k , which become strong within the dislocation core. The component of Ao of 

emergent electromagnetic field is not forbidden by any symmetry. Therefore, it appears and gives rise to emergent 
electric potential (either attractive or repulsive) within the dislocation core. 

Notice, that the simple model of Weyl semimetal with cubic symmetry has been considered in [35] . The Dirac 
semimetal (with cubic symmetry) may, in principle, be described by the two copies of the model of [35]. 
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III. EMERGENT MAGNETIC FLUX CARRIED BY THE DISLOCATION 


In order to calculate the emergent magnetic field we should use integral equation 

y£ijk [ tPdx^ A dx k = f A k dx k , 

2 Js JdS 


( 21 ) 


where the integration goes over a surface in the transverse plane which includes the position of the dislocation. For 
the considered solution of elasticity equations (fl7|) we represent the right hand side of this expression as follows 

( 22 ) 


f A k dx k = VKi 0) + -P ijk f u^dx 1 
JdS a JdS 


IdS a JdS 

The first term in this expression gives the following singular contribution to magnetic field: 


#s fc ingO) ~ bl K- 0) [ dy k (s)6 < ' 3 \x - y(s)), 
ji° 


(23) 


where the integration over y goes along the dislocation axis 1°. 

One can check that the solutions of elasticity equations give ui k ~ 1/r at r —► oo. Therefore, the integral along the 
circle C r = 8S at r — > oo (with the dislocation at its center) in the second term of Eq. (22) gives finite contribution 
to the normalized total flux of the singular gauge field H S i ng : 


where 


iw = fr = 17 ^ 


$n = 27T 


(24) 


(25) 


is the elementary flux (in out units the electric charge is unity e = 1). At the same time the function $(oo) — <3E>(r) 
takes its maximum at r = 0 and decreases fast out of the core of the dislocation. 

In the considered crystals there exist several exceptional vectors Gi (i = 0,1, 2,...), which generate the symmetry 
of Brillouin zone, i.e. momenta k and k + Gi are equivalent. The unperturbed Fermi point is directed along Go and 
is also defined up to the transformations K ^ —> K ^ + Gi. This corresponds to the change of the magnetic flux by 

A<& = b ■ Gi = 2t tN , TV e Z . (26) 


Such a change of the magnetic flux is unobservable for Weyl fermions and Eq. (261 is thus posing certain restrictions 


on the choice of the Burgers vectors. For example, for the layered hexagonal structure of Na and Bi atoms in the 
compound NagBi we have 


4t r „ 
G i = —x, 
3a 


4tt /1 >/3 „ 

G3= u(2 X + ^ V 


4tt (1 „ ^3 

= - — cr. — - 




y 


da 


(27) 


Here a is the interatom distance within each layer in the plane orthogonal to Go || and the material parameter 
£ determines the interlayer distance a z = ||. Due to the hexagonal (honeycomb) structure of the NagBi layers in 
xy plane, we may construct the Burgers vectors similarly to the case of graphene 34J which has also the hexagonal 
structure. Condition (26) gives us the following general expression for the Burgers vectors: 


b = ^2 Ni-rrii 


(28) 


where Ni € Z. The vectors m.j 


3a, 

mo= 2C Z ’ 

, , 3a ^ 73 a „ 

mi = -h + l 2 = ~2~ x + ~^~y i 

m 2 =l 3 -l 2 = —73 ay , 

, 3a ^ 73 a 

m :i = l- { - = ——x + —^~y j 


( 29 ) 


2 
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are constructed from the unit vectors 1 1, 1 2 and I 3 which correspond to the nearest-neighbor Na-Bi bonds of the 
honeycomb lattice in the transverse planes of NaaBi: 


l\ = -a®, 

h=a Cr +) Y y )' 

, /l* \/3.' 

3 ~ a \ 2 X 2 ^ V 


(30) 


For a screw dislocation perpendicular to the layers of NaaBi the displacement vector is given by Eq. (13). The only 
nonzero components of the corresponding deformation tensor (15) are 


,3a ( 


b 3 e 3ab x b j 


(®_l) = u a 3 (x ± ) = 


(31) 


and the emergent electromagnetic field © is given by the following expression: 


P ..3i 


Ai = — V AuK) -\—u 

a 


F 3j 

—e 3ij u J 


An — 0 , 


(32) 


differs from that of Ref. [55] , 
of the emergent magnetic field 

H: 


with some material-dependent constants /3 and (3 Notice that our expression (32) 
Equation (32) leads to the following expression for the (normalized) magnetic flux (24 


- . . K (0) b $ 

$(°°) = — 


27r 


47 r a 


(33) 


For example, in NaaBi the value of K ^ ss 0.26— z, where a z is the lattice spacing in z direction JJ [46]. The value 
of &3 = Na z is proportional to a z . Therefore, the topological contribution to magnetic flux of the dislocation is 
KLAb ~ °' 2 2 ^ N ~ 0.131V. Following an analogy to graphene, where Gruneisen parameter /? ~ 2 we may roughly 
estimate the second term in Eq. (33) as ~ 0.2 N. Then the emergent magnetic flux incident at the dislocation, 
presumably, reaches the value of 27 t at IV ~ 30. 

As it was mentioned above, we may neglect the zero component of the emergent electromagnetic field Aq at large 
distances r>a, where the elasticity theory works. However, such a potential may be present within the dislocation 
core because of the essential change in the microphysics at the interatomic scales. Thus we assume the existence of 
either repulsive or attractive potential 


A 0 { x±) = Vfv 1 /3 cf)(x ± ), 


(34) 


localized at the dislocation core. We neglect possible appearance of such potential far away from the dislocation axis. 


IV. FERMION ZERO MODES PROPAGATING ALONG THE DISLOCATION 


Momentum of quasiparticles that are described by the action of Eq. Q should be much smaller than 1/a. At 
the same time the emergent gauge field A within the dislocation core may be as large as ~ 1/a. Therefore, the 
contribution of emergent gravity to Eq. Q is always small compared to the contribution of the emergent gauge field. 
Nevertheless, for the completeness we consider this contribution in Appendix A. 

In this section we show that the dislocations in a Dirac semimetal host a topologically protected massless 
(quasi)fermion mode, which propagates along the dislocation with the Fermi velocity. This fermionic mode cor¬ 
responds to the zero mode of the transverse Hamiltonian (103). The appearance of the propagating mode localized 


in the vicinity of the dislocation is also known to emerge in topological insulators with lattice dislocations |45| . 

Let us neglect the emergent gravity due to its weakness and concentrate first on the case of the screw dislocation, 
when A 3 = 0. We may apply the gauge transformation, which brings the gauge field to the form 


Ai — Cap' b)j f (x ±_), 

where / is a function of transverse coordinates. Then the Hamiltonian ( [98]) receives the form: 

H (R) = VpV 2/3 a 3 h 

+ VpV i/ 3 n (R) 


(35) 


(36) 
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with 


^ V a (pa+A a (x)y 


a—1,2 


The zero modes of the transverse Hamiltonian (37) are defined as the solutions of equation 

H^ip = 0. 


(37) 


(38) 


Xi=rcos9 , X 2 = rsin$, 
i „ 


p r = —id r 


P9 = ~~ 


Equation (38) is well known in particle physics as it determines zero eigenmodes of a fermion field in a background 
of an abelian vortex m • The magnetic flux of the abelian vortex is equal to the quantized vorticity number n. 
There are exactly \n\ isolated, linearly-independent, zero-energy bound states. These bound states are topologically 
protected by index theorems. For the sake of convenience here we repeat below the derivation of Ref. |TT1. 

We represent ip = e~ a ^ip and rewrite the Hamiltonian in the polar coordinates r, 9 in the transverse plane of 
x±_ = (x\, x 2 ) using 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


and the radial sigma matrices: 


~ m — 

Then equation for the function ip is Tl\ ip = 0, where 


( 0 e~ ie \ 

n i 

( 0 

- ie~ iB \ 

\e ie 0 ) 

q 

II 

( ie ie 

0 ) 


n 


( R ) 


0(r, 0) + a r p r + <T e p e , 


or 


with 


n 


( R) 


e ie n (R) 


+ 


cp{r,9) 

o —i 0 n/(R) 


0(r, 9) e l6/ H: 

'H [ ± ) ft! Pr ± ipe 




n 


(R) 


In the absence of the electric potential <p(r , 9) the zero modes (if they exist) have a definite value of the spin 
projection s = ±1/2 on the z axis. At large r the corresponding coordinate parts of their wave functions satisfy the 
relations 


(Pr ± *Pe)V’± ° = 0 • 


f r - dr 

/M)= / Hr, 9)-, 

./n T 


Next, we chose 


so that the angular d-component of the gauge potential (|35j) gets the following form: 

, Hr, 9) 

Ap = -. 


(45) 

(46) 

(47) 


The axial symmetry of the problem implies that at large distances r the function >I>(r, 9) is independent of the polar 
angle 9. Therefore, at large r the function $(r) is the magnetic flux within the circle S r of radius r: 


$(r) = 




1 


Cijkdx 3 A dx k H\x , y ), 


(48) 


where the surface S r belongs to the plane which is orthogonal to the dislocation. We come to the following solutions 
of Eq. (381 for the right-handed zero modes [J3] : 


ip ( ±\r, 9) ~ r m e ±imB ^Jo ®( r ’ e ') d r 


( 49 ) 
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where the integer m is the angular quantum number. 

The solutions (49) are localized in a small vicinity of the dislocation core provided that the angular quantum number 
satisfies the following condition 


m — 2s<f> < —1, 


(50) 


where 


» <J> 

$ = — = lim d>(r, 9 ), 

<f>0 r—^oo 


(51) 


is the total flux of the intrinsic magnetic field H normalized by the elementary magnetic flux (25). If Eq. (50) is 
satisfied then the corresponding probability distribution is convergent at large r: 



rdrd9\ip\ 2 


= 1 . 


(52) 


Notice that the ultraviolet cutoff £ is of the order of the lattice constant a. 

In addition, there exist two solutions of Eq. (45), which may not be normalized and which have their maxima at 
the dislocation core provided that 


TO = [2s6(oo)], 


(53) 


where [2s<f>] is the integer part of 2s4>, which is the maximal integer number that is not larger than 2s<3E>. 

The probability distributions of the considered solutions are convergent at small r for m > 0. Therefore, in the 
absence of both the nontrivial vielbein and the electric potential, the zero modes that are not singular at r —> 0 and 
are not localized on the boundaries of the system, should satisfy 0 < to < 2s$(oo). Such modes exist for s<h(oo) > 0 
and are enumerated by the values of orbital momentum 

ro = 0,...,[2s$], (54) 


We neglected in this derivation the potential <f>. However, it is localized at the dislocation. Therefore, the zero modes 
in the presence of electric potential (if they exist) have the form of Eq. ( |49| ) at r > o. Recall, that Eq. Q works 
for the momenta of quasiparticles much smaller than 1/a. Therefore, the solutions of Eq. (38) localized at the 


dislocations, presumably, do not represent physical zero modes. The only solution that remains is the one with 


to = [2 s<l], s = -sign<I> 


(55) 


Fortunately, the field (f> cannot affect the energy of this solution because the probability density corresponding to this 
solution of Eq. ( [45] ) is dominated by the distances far from the dislocation core, so that we can neglect completely the 
region of the dislocation core. The vielbein for this solution also gives small corrections compared to the contribution 
of emergent magnetic field. Therefore, the strong gravity and the potential </> at r ~ £ cannot affect the main properties 
of this solution: it certainly survives as the zero mode and still has the definite value of the projection of spin to the 
z axis. 

In the case of edge or mixed dislocation we should take into account the appearance of a nonzero third component 
of the emergent gauge field: 


A 3 (x±) 


~2 (j3ib±xx 


■ d 2 .r 


)■ 


Then far from the dislocation core one gets 


H 



e ie nf ] vA 3 (r,6)\ 
-iyA 3 {r,0) j 


(56) 


(57) 


Perturbative corrections to the eigenenergy due to the presence of A 3 may be nonzero, in principle, but for the 
mode with to = [2 s<f>(oo)] those corrections may be neglected because all integrals are dominated by the regions with 
r —>■ oo while A 3 ~ 1/r at large distances. 

Thus we come to the conclusion, that the only zero mode existing around the dislocation is the one with 

s = -sign4>(oo). 


to = [2 s$(oo)], 


(58) 
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(R) 

The zero mode (49), (58) of the transverse Hamiltonian Ti _j_ ' corresponds to the zero mode of the full Hamiltonian 
provided the longitudinal momentum is zero p 3 = 0. At the same time it corresponds to a linear branch of 
spectrum of the full Hamiltonian Ti ^ with the corresponding dispersion law: 


E {R) » u FF 2 / 3 sign(4>)p 3 . 


(59) 


This branch crosses zero energy level at p 3 = 0. 

Similar considerations can also be applied to the left-handed Hamiltonian / HS l \ where the only physical zero mode 
of the corresponding transverse part T~L^ is 


(r, 9) ~ r m e i2sm9 - 2s So Hr) w 


(60) 


with the quantum numbers 


m = [2s$], s=-sign<3E>. 

This mode corresponds to the branch of spectrum with the dispersion 

ss —w F zz 2 / 3 sign(<i>)p 3 . 

The right-handed and left-handed fermionic modes propagate along the dislocation with the velocity 

vr = —v L = u F i' 2/3 sign(<l), 


(61) 


(62) 


(63) 


which is nothing but the corresponding component of the anisotropic Fermi velocity. Thus, the right-handed massless 
quasiparticle propagates up or down along the dislocation depending on the sign of the flux >!>. The left-handed mode 
always propagates in the opposite direction compared to the right-handed mode. 

Sure, these mode may be coupled with the counter-propagating zero modes, for example, by the intra - valley 
and/or inter - valley impurity scattering. This effect is neglected here. But this is what is also neglected always in the 
Dirac semimetals, when one speaks of the massless Dirac fermions emergent in those materials. From the theoretical 
side the gapless nature of the emergent Dirac fermions in Dirac semimetals should be protected by momentum space 
topology: the symmetry - protected topological invariant in momentum space should exist that protects the Dirac 
cones (see [36]). The corresponding symmetry is specific for the crystals of the semimetals. For example, for the case 
of Na 3 Bi the effective model (that follows from the microscopic electronic model) is given in [46]). This model utilizes 
the time reversal, inversion, and D^ h symmetries. Their composition (possibly, supplemented by something) should, 
in principle, protect the Dirac cone for the case of Na 3 Bi. But the extensive theoretical study of this issue is out of 
the scope of the present paper. 

Anyway, the symmetry responsible for the gapless nature of the modes localized around the dislocation is the same, 
which is responsible for the gapless nature of the bulk quasiparticles. This symmetry is not exact, and it is indeed 
broken softly by the impurity scattering and also by the collisions with the change of chirality. The corresponding 
parameter is the mean free path, which serves here as the infrared cutoff of the theory and as the typical size of the 
region, where the zero modes of the present section are localized (see discussion in the forthcoming sections). 

Notice that Eqs. ( |59| ) and ( |62| were derived in the assumption that the magnetic fluxes of the emergent magnetic 
field for the right-handed and the left-handed $ F quasiparticles are the same, = <£>. However, if in a 

Dirac semimetal the constants /3, ?F differ for the left-handed and the right-handed fermions, then the corresponding 
gauge fields (17) are also different, and in this case the magnetic flux entering Eq. (59) will be different from the flux 
in Eq. (62). 


V. CHIRAL ANOMALY IN DIRAC SEMIMETALS ALONG THE DISLOCATION 

In the presence of external electric field E the states that correspond to the described above zero modes flow in the 
correspondence with the following equation: 


<P3> = E 3 . (64) 

Now let us take into account, that the studied model has the infrared cutoff l/£o, where £o a. Then the 
zero modes of TT^ and TL^\ and the corresponding branches of spectrum of the propagating modes of the full 
Hamiltonians Ti^ and Ti^ obey the following properties: 
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1. The propagating fermion modes are not localized at the dislocation core. Instead, the region of space around 
the dislocation of size £o dominates, where l/£o is the infrared cutoff of the theory. 

2. The propagating fermion modes have the definite value of the spin projection on the dislocation axis: s = ^sign $. 
The corresponding branch of spectrum for the right- and the left-handed fermions is given, respectively, by the 
following dispersion relations: 


£r/l{P3 ) = ±2sv f v 2/3 p 3 . 


(65) 


3. The propagating mode appears for any dislocations including 

than nnit.v 


than unity. 

The total production of the right-handed quasiparticles per unit length of the dislocation is given by: 

En - 

Qr = ^rSign$, 


( 66 ) 


where the unit vector n is directed along the dislocation. In the following we assume for simplicity, that the signs of 
the emergent fluxes $ experienced by the right-handed and the left-handed fermions coincide in the Dirac semimetal. 
Therefore, the production of the left-handed quasiparticles in Dirac semimetal is given by 


Ql = ~qR ■ 

The production of the quasiparticles may be written as the anomaly in their currents 

3l = I e 0*01 j l( x ) » Jr = I e(x)\J^{x ), 


J — 3 R + 3 L ■. 


3 5 3 R 3 L 


(67) 


( 68 ) 

(69) 


where the covariant currents are defined according to Eqs. 0 and ([8]). In a local form the anomaly may be expressed 
as follows: 


FjYI a 

(d„j£(x)) = - fo{x±) sign $ . 

7r 


(70) 


where the function f 0 can be read from Eqs. (49), (60): 


fo(x±) = 


exp (—|*l|/£o) (f^) 


—2(|<S[-[|4>|]) 


27r^r(-2|$|+2[|$|]+2) 


(71) 


This function is normalized in such a way, that 


27t / rdr/o(r) = 1. 


(72) 


The quantity l/£o has the meaning of the infrared cutoff of the theory and the factor exp (— |:e_l|/£o) appears as the 
infrared regulator. We imply that the size of the semimetal sample is much larger than the infrared cutoff £oj while 
£o is much larger than the size of the dislocation core £ ~ a, £o £■ Thus the chiral anomaly due to the zero mode 
with m = [|$|] is localized within the tube of size £o centered at the dislocation. 

In our derivation we follow the procedure first applied by Nielsen and Ninomiya in [ 48] . This derivation is now 
standard and it was repeated many times. In particular, the derivation of Eq. (66), (67) may be found in Sect. 18.1.1 
- 18.1.3. of [86] . Eq. (70) is the local form of Eqs. (66), (|67|. It appears as the local violation of the continuity 


equation that follows from the appearance from vacuum of the state with the wave function of Eq. (601. 


It is worth mentioning, that for a single dislocation, if for some reasons the contributions to the emergent magnetic 
flux of the dislocation due to (3ijk in E q. (fTf] may be neglected, the typical values of the Burgers vector are such that 
|$| < 1. Therefore, according to Eqs. (|58faiid (61), for a single dislocation the zero mode corresponds to the angular 
momentum m = 0. 

We may rewrite the expression for a chiral anomaly caused by a single dislocation in a Dirac semimetal as follows 


(d/Js (*)> = 


EB 

2 n 2 


(73) 
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where the effective magnetic field B responsible for the chiral anomaly is given by 

B(x±) = 27rn/ 0 (a;_ L )sign4>, (74) 

where the function /o is given by Eq. Sl¬ 
it is worth mentioning, that the above consideration refers only to those branches of spectrum, which are described 
by the low energy effective field theory. In the presence of electric field the pumping of the quasiparticles from vacuum 
may also occur at another branches of spectrum. Ideally, this pumping process should be treated with the help of a 
microscopic theory and is out of the scope of the present paper. 


VI. THE CASE OF DISLOCATION DIRECTED ARBITRARILY 

In this section we consider the dislocation directed arbitrarily. Without loss of generality we consider the dislocation 
directed along an axis, which belongs to the (yz) plane. The angle between the dislocation and the z axis is denoted 
by ip. Let us rotate the reference frame in such a way, that the z axis is directed along the dislocation. In the new 
reference frame the tensor / has the form: 

/ is- 1 ' 3 0 ° \ 

/ = J 0 v~ 3 ! 3 cos ip i/ _1 / 3 sin tp J (75) 

\ 0 —i/ 2 / 3 sin ip v 2 / 3 cosip J 

Let us apply the transformation of spinors i/j —► e l ^ a with tga = vtgip. In the transformed frame the tensor / is 
modified: 


f = 


( v - x ' 3 0 

0 i/ 1//3 \/j/ -4 / 3 cos 2 i p + i/ 2 / 3 sin 2 ip 


V 


0 

(1 — v 2 ) sin 2 ip 


\ 


Ivy/ v~ 4 / 3 cos 2 cp-\-i/ 2 / 3 sin 2 ip 


\Jl/~ 4 / 3 cos 2 <p+i/ 2 / 3 sin 2 <p J 

The one-particle Hamiltonian for the right-handed fermions becomes as follows 

n (R) = v F f 3 a 3 p3 + v F fi<j 2 p 3 + vpfl'H^ 

with 


n 


(R) 


^{Pl - Ai) + 4j-cr 2 (p2 - A 2 ) 

Jl 


f 3 

~J[ a 3A3(x_ l) + (j){x_ l) 

Jl 


Now we perform the coordinate transformation 


y 


fj 

fl 


y , 


Ay "A ~ 2 A y , 
J 2 


(76) 


(77) 


(78) 


(79) 


and notice that the equation for the zero mode of the Hamiltonian T-L± becomes the same as the one discussed in 
Section IV Thus we arrive at the expression for the anomaly in quasiparticle current of Eq. (66). The resulting 
expression for the chiral anomaly in the Dirac semimetal is again given by Eq. (73). 


VII. CHIRAL DENSITY AND CHIRAL CHEMICAL POTENTIAL AROUND THE DISLOCATION IN 

THE PRESENCE OF ELECTRIC FIELD 


The Dirac senrinretal possesses two cones, each of which hosts one right-handed and one left-handed Weyl fermion. 
Since the processes operating in these two cones are equivalent, we concentrate on one cone hereafter taking into 
account the fact of the double degeneracy later. 

The evolution of the local chiral density around the dislocation is governed by (i) the generation of the local chiral 
charge due to quantum anomaly at the dislocation given by Eqs. (73) and (74), (ii) the spatial diffusion of the chiral 
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charge and (iii) the dissipation of the chiral charge density ( [85] ) . In the other words, in the presence of the external 
electric held, the zero modes, distributed around the dislocation and propagating along the dislocation, accumulate 
the chiral charge. The accumulated chiral charge diffuses (due to, basically, thermal diffusion and scattering) and also 
dissipates (due to the chirality-changing processes) around the dislocation. We estimate these effects below. 

At the distances larger than the size of the dislocation core r > ( ~ a we may neglect the presence of the emergent 
magnetic held H of the dislocation. Therefore, in order to relate the chiral chemical potential p$ with the chiral 
density p 5 at finite temperature T we use an approximation, in which the relevant modes of the quasiparticles are the 
plane waves of the continuous spectrum. Thus, we neglect gauge held completely and calculate the thermodynamical 
potential (see [3] and references therein): 


fi = r E E 

S=±l c=±l 


d 3 p 

(2tt) 3 



W P,S+C (45 \ 

T )’ 


(80) 


where the quantity c = ±1 labels right- and left-handed chiralities, s= ±1 is the projection of spin (multiplied by 
two) to the auxiliary vector k a (p) = faPi, while pi is momentum of the quasiparticle. The chiral chemical potential 
p, 5 is the difference between the chemical potentials associated with the fermions of right-handed and left-handed 
chiralities: 


d5 = 2 (dR ~ Pl) ■ 


In Eq. (801 the dispersion of the quasiparticles in terms of the vectors p and k is as follows: 


(81) 


Up, s =CSV F J fafaPiPj = CSV F \k(p)\, 


(82) 


where v F is the Fermi velocity (63) which enters the dispersion relation for the chiral fermions (65). The matrix / is 


given in Eq. (12). For the momentum parallel to the z axis, p = (0,0,p 3 ), Eq. (82) leads to the dispersion (65). In 
our calculation we work (following, e.g., Ref. [3]) in the adiabatic approximation assuming that the chiral cl 
potential is a slowly varying function of space and time. 


remical 


Since the determinant of the matrix / is equal to unity, we get for the thermodynamical potential (80): 


n = t 


E E 

s=±l c=± 1• 


d 3 k 

(2 


log (l 


I k | +C/X 5 


(83) 


The chiral density is given by the derivative of the thermodynamical potential with respect to chiral chemical 
potential /i 5 : 


till 1 x -> x -> 

dp 5 r> ' n ' 2 2— 1 2— 1 


ck 2 dk 


v F k + cp .5 


(84) 


s=±l C—it 1 ' 


Obviously, for p 5 = 0 the chiral density vanishes. For \p 5 \ <C T we may evaluate the term in p 5 linear in p 5 
differentiating the last expression with respect to p 5 : 


P5 


P5 

2 tt 2 T 

2 p 5 


E 


~k 2 dk 


S=±1,C=±1 ° 

[ k 2 dk 


(i+e-^y 

PsT’ 


7 t 2 T J cosh 2 (2^#) 3u| 


(85) 


Equation (85) is valid provided that certain conditions are satisfied. First, the size of the dislocation core £ ~ a 
(where the low-energy physics, and, consequently, the thermodynamic relation (85) both become inapplicable), should 
be small compared to the wavelength of the typical thermal momentum At ~ 1 /p F ~ v F /T that contributes to 
relation (851. According to our estimates (see below, Section VIII) this condition is satisfied even for the room 
temperature T ~ 300K (with corresponding Ay ~ 4 x 10 -8 m) because £ ~ 10” 9 m <C At according to Eq. ( 114 ). 
Second, the magnetic field of the dislocation should not affect considerably the plane waves contributing to Eq. (85). 
To this end one can consider a wavefunction of a particle that circumferences the dislocation, and compare the 
contributions to its phase coming from the magnetic field and from the usual kinetic factor e lpx . The former reaches 
its maximum at r = £ 0 , being equal to the total flux = $ ~ 27 t while the later can be estimated as 54>t = ‘^'RPtS.o- 
Thus, the second condition requires 8 <j>& <C 54>t or At <C Co which is also satisfied according to Eq. (113). 
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The nonconservation of the axial charge can conveniently be written in the following form (see [3]): 


dp5 

dt 


+ Vj 5 = 


P 5 
T V 


1 

27T 2 


BE, 


( 86 ) 


where the first term in the right hand side corresponds to the dissipation of the chiral charge with the rate given by 
the chirality-changing scattering time ry while the second term describes the generation of the chiral charge due to 
the quantum anomaly around the dislocation. In fact, this equation appears also in the chiral kinetic theory developed 
relatively (see, for example, |491 150] 1 when the kinetic equations are taken in the relaxation time approximation (also 
known as r approximation). 

The chiral current, 


j 5 = ~D 5 Vp 5 . 


(87) 


is given by the diffusion of the chiral charge p 5 with the corresponding diffusion constant D 5 . This equation implies 
that the chiral charge is able to diffuse. We do not have the microscopic description of the diffusive nature of the 
chiral charge, and at the present level it is considered here as the hypothesis. 

We assume that the Dirac semimetal has zero usual chemical potential for the Dirac quasiparticles, p = 0. Moreover, 
we consider a linear approximation so that the transport effects, which are discussed here, do not generate a nonzero 
p. Substituting Eq. (f87l) into Eq. (| 86 | one gets the following equation for the chiral charge density: 




( 88 ) 


In the constant electric field, dE/dt = 0, the chiral charge P 5 relaxes towards equilibrium dp 5 /dt = 0 at late times 
t t y. The equilibrium chiral charge density is given by a solution of Eq. ( 88 ) with vanishing left hand side: 


P5(x) = 


1 


2tt 2 T>, 


d 3 yG^(x~y,X){B(y)-E(y)), 


where 


(-A + L v 2 ) G (3) (x -y; A) = S(x - y ), 
is the three-dimensional Green’s function and 

Lv = VD5TV , 


(89) 

(90) 

(91) 


is a characteristic length which controls spatial diffusion of the chiral charge. 

Working in linear approximation we consider weak electric field E, so that the chiral imbalance can always be 
treated as a small quantity, /x 5 <C T, so that the linear approximation in Eq. ( |85[ ) is justified. In the absence of the 
usual chemical potential p, one gets from Eq. (85) the following relation between the chemical potential and the chiral 
charge density (89): 


/ \ 3 V rp . . 

M5W = -J^-Pi-o{x) . 


(92) 


Thus we see, that the dislocation produces the chiral charge which spreads around the dislocation, effectively creating 
an excess of the chiral chemical potential at the characteristic distance Lv from the dislocation axis. Notice, that 
this chiral chemical potential corresponds to a single Dirac point with a pair of Weyl fermions. 

In the above derivation we neglect the gradient of temperature. This may be done for sufficiently small external 
electric field, when the temperature remains almost constant at the characteristic length of the problem that is £ 0 - 
Very roughly, in equilibrium the heat generated by electric held ~ crE 2 (where a is the total conductivity that includes 
Ohmic contribution) should be equal to the divergence of the heat how kVT (where k is the thermal conductivity). 
For our estimate we use the Wiedemann-Franz law k ~ aT. This gives for the characteristic length (at which the 
temperature is changed considerably): 


1 AT E 2 
& ~ ~T~ ~ T 2 


(93) 


According to our estimates (see below Sect. VIII) at room temperatures the condition £0 leads to \E\ <C 1 

V/cm. This condition provides, that temperature remains constant within the region of size £0 around the dislocation. 
However, temperature may vary within the whole semimetal sample if its size is much larger than £o- 











16 


Now let is consider practically interesting case when the dislocation is a strait line centred at the origin, X\ = X 2 = 0 
and directed along the X 3 axis. The dislocation induces the intrinsic magnetic field 


B(x±) = B z (x±)n 1 


(94) 


which is also directed along the X 3 = z axis (here n is the unit vector in z direction). The intrinsic magnetic field is 
a function of the transverse coordinates Xx = (^ 1 ,^ 2 ) which takes nonzero values in a (small) core of the dislocation. 
In our model approach we consider the field given by Eqs. (71) and (74) 


B z ( x±) = sign 4> 


exp (—|x l |/£o) (f^) 


-2(|$[-[|$ID 


where 


e 0 2 T(-2|d>|+2[|4>|] + 2) 


a=|$|-[|$r 


(95) 


(96) 


is the fractional charge of the normalized flux <1> |. The effective magnetic field (95) is distributed around the dislocation 


with the characteristic length £0 that is much larger than the size £ of the dislocation core (the latter is of the order 
of a few lattice spacings a). In Fig. ^ we show the field (951 for a few values of the fractional part of the flux a. 



FIG. 3. The effective intrinsic magnetic field B z (x j_) in Eq. (|95| vs. the distance from the dislocation core x± plotted for a 


few values of the fractional part of the absolute value of the normalized flux (511 a = 0.1, 0.5, 0.9. 


This effective magnetic field carries the unit of the elementary flux (25): f d 2 x±B(x±) = 27rsign<f>. As we have 
discussed earlier, this effective magnetic field is associated with the propagating zero modes bounded at the dislocation. 
In Eq. (95) the total flux *I> of the intrinsic magnetic field H is of a geometrical origin. The flux is a quantity of the 
order of unity (51), in terms of the elementary magnetic flux ((251). For the straight dislocation (95) the axial anomaly 


generates the axial charge which spreads in the semimetal in the transverse directions according to the equilibrium 
formula (l89|). 


For the completeness we represent the numerical estimates for the encountered above constants that characterize 
the semimetal in Appendix B. 


VIII. CONCLUSIONS 

In this work we discussed certain effects of anomalies in the Dirac semimetals NaaBi and CdaAs 2 caused by the 
dislocations in their atomic lattices. This chiral anomaly is operational without any external magnetic field unlike the 
conventional chiral anomaly that was discussed for Dirac semimetals, for example, in [4] . The dislocation appears as 
a source of the emergent magnetic field as it carries the emergent magnetic flux. The emergent flux gives rise to the 
zero mode of the one-particle Hamiltonian for the fermionic quasiparticles. The fermionic mode is a gapless excitation 
which propagates along the dislocation being localized in the area of the size £0 around the dislocation (here l/£o is 
the infrared cutoff of the field theoretical approximation used in our approach). The length may also be identified 
with the mean free path of the quasiparticles. For example, for Cd 3 AS 2 it is of the order of 200 gm. This propagating 
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zero mode corresponds to the branch of spectrum of the quasiparticles with the spin directed along the magnetic flux 
of the dislocation. 

In the presence of an external electric held the spectral how along the zero-mode branch of the spectrum leads 
to the pumping of the quasiparticles from the vacuum. Since the right-handed and left-handed quasiparticles are 
produced with opposite rates, the pumping process corresponds to the chiral anomaly (73). The production rate of 


the chiral density is controlled by a scalar product of the usual (external) electric held E and the effective (internal) 
magnetic held B carried by the dislocation. One should stress the following subtle fact: the dislocation carries also 
the emergent magnetic held H [given in Eq. (33) for the example case of NaaBi] which gives rise to the existence 
of the mentioned propagating zero mode. However, the emergent held H does not directly contribute to the chiral 
anomaly (73): it is the effective magnetic held B that is expressed via the density of the zero mode - which 
enters the anomaly relation (73). 


In the other words, the emergent magnetic held H with the hux (33) leads to appearance of the right-handed (49) 


and the left-handed modes (60) localized at the dislocation and propagating along its axis. The external electric held 
E, parallel to the dislocation axis, produces the chiral charge by pumping these modes form the vacuum at unequal 
(in fact, opposite) rates proportional to the scalar product EB. The process can be formulated via the chiral anomaly 
equation (73]), in which, however, the magnetic held B is expressed through the density of the wave functions of the 
mentioned zero modes (with B ^ H). 

The chiral anomaly gives rise to a nonzero charge density localized around the dislocation axis with the characteristic 
localization length Ly ~ 2 /iin for Cd 3 As 2 - The slowly varying chiral density can be expressed via a (space-dependent) 
chemical potential. 

In principle, there are various ways to create a Dirac semimetal with dislocations. In general, the growth of an 
atomic crystal may be organized in such a way, that the dislocations appear along the chosen direction with the 
chosen values of the Burgers vector. In addition, dislocations may also appear as a result of plastic deformations of 
the crystals m- This opens a possibility to observe the effects of the dislocation-induced anomaly experimentally. 
The chiral density that is formed around the dislocation in the presence of external electric field should affect transport 
properties of the semimetal. In order to calculate the corresponding observable quantities, this is necessary to use 
kinetic theory modified accordingly in order to take into account the appearance of the chiral density around the 
dislocations driven by chiral anomaly. However, the solution of this problem is out of the scope of the present paper. 
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APPENDIX A. EMERGENT GRAVITY AROUND THE DISLOCATION 


In this section we briefly consider emergent gravity around the dislocation. Let us represent the action for the 
right-handed fermion in the following way: 

Sr = j d 3 xdt^f{x,t) \e(x)\e^(x)id t 

-\e(x)\e k a (x)a a o (P k - A k ) ^(x,t) 


= / d 3 x dt^(x,t) 


Wt - n {R) 


T(x,t), 


where P k = — i\7 k , = y / \e(x)\e 3 (x) '3/, a = 0,1, 2, 3 and k = 1, 2, 3. The one-particle Hamiltonian is given by 


where 


= f k (x)a a o (P k - A k ), 

e a( X ) 


/£(*) = 


) 0 ) ' 


(97) 


(98) 


(99) 


We used here the following chain of relations: 
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d 3 x dt 


d 3 x dt 


^ '(x,t) 


V\ e ( x )\ e o( x ) 

1 


W{x)\e k (x)a a di 




V\ e ( x )\ e o( x ) 




V\ e ( x )\ e o( x ) 


^{x,t) 


e{x)\e° Q {x) 


\e(x)\e^(x)a a di^(x,t) - d^(x,t) |e(x)| e k (x)& 


e(x)\e k a (x)a a 

1 


^(x,t) 

V\ e ( x )\ e o( x ) 

V{x,t) 


e(x)\e° 0 (x) 

x*dt {nx,t)ft(x)^Mx,t) - [«*,*)] s 2/ XxdtiM Kw^.»] $(*,*)■ 


( 100 ) 


We represent /£(#) as follows 


fa (x) * f k - rtde b a (x) 


fo( x ) ~ -VFfbSe b 0 {x), a,b,k = 1,2,3, 


( 101 ) 


where the expressions for the small variations of the vierbein held <5e^ can be read off from Eq. (16 1 . 


by 


The one-particle Hamiltonian for the right-handed fermions in the presence of a dislocation along the z axis is given 


= v F v 2 ' 3 a 3 p 3 - v F v 2 ' 3 X 

+v F v~ 1/3 nf ) , 


a —0 


( 102 ) 


where the transverse part of the Hamiltonian is 


Hi® ~ Y “ X v a 5 e a( x ±) °Pk 


/c—1,2 


+</>(*, y) - w 3 ^3(®±) 


(103) 


[o- 3 fc 3 (x_L) +Seo(x ± ) 


°Pk ■ 


k= 1,2 


In a general form, the dislocation-induced deformations of the vielbein held 6 e% in the Hamilt onian (102), (103) 
can be expressed via components tensor 7 of Eq. (16) and the relations given in Eqs. (99) and (101). However, in 
certain symmetric cases the form of the Hamiltonian may be simplihed. Consider, for example, the case, when the 
screw dislocation is directed along the 0 axis of the NasBi atomic lattice (or, equivalently, along the vector K < - 0 '). 
Then, one can write the following expression for the deformations of the vielbein: 


Se k a = lx K ak ^u 3 ^ + -y 2 K akj u 3j , 


Se 3 = 7363 kjW 


31 


■ 74 W 


3k 


k = 73t; 

Sel = 7 5 u 3k + 7663 kjU 3j 

5el = 0 , 

Se k = 0, 


.3k 


a, k = 1, 2, 
k = 1 , 2 , 
k = 1,2, 

k = 1,2. 


(104) 

(105) 

(106) 

(107) 

(108) 


,i3 


Here we have used the fact that the only nonzero components of the tensor of elastic deformations (15) are if 1 = u 
with i = 1,2 given in Eq. (31). Moreover, we took into account that the dislocation is directed along the z axis 
which is perpendicular to layers of honeycomb lattices formed by Na and Bi atoms in the transverse (x, y) plane. The 
requirement to respect the C 3 rotational symmetry of the honeycomb lattice in the (x, y) plane allows us to define 
two tensors from the nearest-neighbor vectors (30): 


J{ ‘jk _ _ 
J^ijk _ _ 


3 a 3 
4 

3 a? 


X vi 1 

6 = 1 , 2,3 

X i ui i 

b= 1,2,3 


b ^3mk 


(109) 

( 110 ) 


which enter Eq. (104) with the material-dependent prefactors 71 and 72 , respectively. The only nonzero elements of 
these tensors are: 


- A ' 111 = K 122 = K 212 = K 221 = 1, 

k 112 = k 121 = k 211 = -k 222 = 1 


( 111 ) 
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The tensor (109) was first introduced in Refs. 


The appearance of the second tensor structure (110) in 


Eq. (104) is a nontrivial fact because the tensor K l: ‘ k is not invariant under P-parity transformation of the 3d space. 
The P-parity odd part is justified, however, by the chiral property of the screw dislocation, because the left-handed 
screws and right-handed screws are not equivalent as they cannot be superimposed on each other with the help of 
rotations only. Therefore, P-parity odd terms may appear in the Hamiltonian. 

Similar arguments lead to the appearance of the other four material-dependent terms in Eqs. (105) and (106) with 


parameters 73,..., 76. Equation ( 108 ) originates from the supposition that the dislocation does not break T invariance 
so that all components of the vielbein involving one tempo ral and one spatial components must be zero. Notice that 
the deformation of the e[j does not enter the Hamiltonian (98) because = 1 according to Eq. (99). 

One can see, that even in this relatively simple case, the expressio ns in Eqs. (|1^4j)-( 108 ) contain six phenomenological 
parameters 7 ,;, and the resulting Hamiltonian given in Eqs. (103) and (|102|), is rather complicated. 


APPENDIX B. CERTAIN NUMERICAL ESTIMATES 


We take for a reference the Dirac semimetal Cd 3 As 2 - The diffusion length of the axial charge for this semimetal 
was experimentally estimated in Ref. [ 7 ] as Ly ~ 2 x 10“ 6 m. This quantity turns out to be almost temperature- 
independent in a wide range of temperatures T = (50 ~ 300) K. A rough estimate of Ref. jS] gives for the relaxation 


time Ty ~ T tr « 2 x 10 10 s. Then from Eq. (91) one finds 


D5 — Ly/ry 


2 x 10“ W/s. 


( 112 ) 


Correspondingly the inverse infrared cutoff £0 may be estimated as follows: 

£0 ~ v f t v ~ ^300 • 10 6 j ■ 2 x 10“ 10 s = 3 • 10 “ 4 m. (113) 

In this estimate we use the value of Vp for Cd 3 As 2 that is around 1/200 speed of light. 

The value of £0 should be compared to the size of the dislocation core 

£ ~ 10 “ 9 m, (114) 

and to the value of the diffusion length 

Ly ~ 2 • 10 “ 6 m. (115) 


Thus we see that in practice the suggested limiting case is indeed realized: 


Ly <C £ 0 


(116) 


and the typical value of parameter x = £o/Ly is x ~ 100. 

Notice, that we used in the present paper the relativistic system of units, in which the only dimensional unit is 
the electron-volt (eV). For example, our distances are measured in eV _1 . We give the estimate in relativistic units, 
where it is expressed through eV or 1/m, where the unit of distance (m) is related to eV “ 4 according to the standard 
relation [200 MeV ] -1 ss 1 fm = 10 “ 15 m. Then the quantities under consideration may be expressed in the SI system 
using the definition of its unit of electric current (A) as Coulomb/s. The SI current equal to one Ampere corresponds 
to the relativistic current equal to 1 /(ec) in the units of 1/m, where e is the charge of electron (in Coulombs) while c 
is the speed of light (in m/s). 

A room temperature corresponds to T ~ 300 K ta 0.025 eV. At the same time D 3 /c ~ 6.7 ■ 10 _11 m = 6.7 • 10 4 fm ss 
3 • 10“ 4 eV _1 . One should also take into account that the typical value of vp in Dirac semimetals is of the order of 
~ 1/200 of the speed of light. We denote by v the degree of anisotropy of the Fermi velocity. In practise in Cd 3 As 2 
[ 2 ] v fJ\ ~ Vpf 2 ~ c/200 while f 3 ~ 0.1/i. In Na 3 Bi [I] Vpj\ « 4.17 x 10 5 m/s, vpf 2 ~ 3.63 x 10 5 to/s ~ c/800, while 
vpf 3 ~ 1.1 x 10 5 m/s. Thus here f 3 ~ 0.27/i. 


Notice, that for small values of the Burgers vector the value of $ may be much smaller than unity. Say, in NasBi 
the minimal topological contribution to $ is of the order of 0.1 (see Section III). However, for larger values of the 
component b 3 of the Burgers vector, the value of $ may always be made close to unity. Besides, the contribution of 
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the second term to the magnetic flux in Eq. (33) also increases the total value of $. The chiral density should also be 
enhanced in a “forest” of dislocations, which are parallel to each other. 
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